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1. INTRODUCTION 

Photovoltaic (PV) based power generation systems are among the key players in the world’s future 
energy mix due to their cost-effectiveness and the advanced ability to connect them to the grid [1]. However, 
some issues still need to be addressed to maximize the output energy from the PV-based power generation 
systems. One of the emerging issues is the extraction of the maximum power under non-uniform weather 
[2]—[4]. The hill-climbing and its related maximum power point tracking (MPPT) algorithms, such as perturb 
and observe (P&O) [5]—[7] and incremental conductance (INC) [8]—[10] are simple and powerful methods in 
locating the maximum power point under uniform conditions. However, under non-uniform weather, they 
may be stuck at a local maximum power point (LMPP) because of the existence of multiple maxima in the 
power-voltage (P—V) curve [11]. 

PV based power generation systems with reconfigurable PV modules were investigated in [12] 
and [13] to minimize energy drop due to the partial shading effect. Although reconfiguration techniques can 
produce higher output power than the fixed configuration in the PV based power generation systems, 
reconfiguring the PV modules’ location is not trivial and needs additional hardware. This will cause higher 
costs and greater system complexity. 
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The artificial bee colony (ABC) is adopted to detect the partial shading conditions (PSCs) and the 
P&O is used to locate the global maximum power point (GMPP) in [14], forming a two-stage MPPT 
algorithm. Another two-stage global maximum power point tracking (GMPPT) algorithm is proposed in [15], 
where particle swarm optimization (PSO) is responsible for detecting the PSC and the P&O is accountable 
for the final GMPP tracking. Statistical evaluations of different MPPT algorithms based on stochastic 
optimization techniques regarding stability, success rate, and convergence speed are done in [16] and [17]. 
However, as in any stochastic optimization technique locating the global maxima is not guaranteed, which 
may cause energy loss. Besides, the time to track the GMPP in these techniques is relatively long. 

Other GMPPT algorithms divide the P—V curve into several voltage segments. Then at each voltage 
segment, the LMPP is tracked using P&O or related MPPT methods. The GMPP is then located comparing 
all LMPPs [18]-[20]. The module open-circuit voltage (V.¥°") of the PV array is usually used for voltage 
segmentation where 0.8 VMOP is adopted in [18] while 0.5 VOP is preferable in [19]. How much fraction 
from the V7 suffices to track the GMPP is unclear, and the actual GMPP can be missed if the voltage step 
is not adequately selected [20]. Also, the transient voltage of the PV generation system must be wisely taken 
care of while scanning the P—V curve; otherwise, misleading LMPPs will be registered instead of the actual 
LMPPs. 

Ghasemi et al. [21] replaced the voltage step scanning with ramp function scanning to mitigate the 
PV voltage transient effect on the MPP. The work in [22] uses the current variation to detect the MPP at 
sampling multiple of VO? , then the hill-climbing algorithm is called when a peak is detected to locate the 
exact LMPP. By estimating an upper power limit for each sub-region at some samples of the I-V curve, the 
search region is minimized in [23]. The procedure will continue until the voltage segment length is less than 
the VOP , then the GMPPT is located using P&O. This method may fail if the peak points are closed to each 
other. 

The algorithms above cannot deal with the hidden points because they rely on power variation to 
detect partial shading (PS). The hidden points are points in the P—V curve caused by a transition between the 
GMPP and LMPP when the PSCs are changed, which cannot be observed by monitoring only the power 
difference of the PV-based power generation system. Therefore, this paper proposes an effective and fast 
algorithm to extract the maximum power from the PV based power generation systems under non-uniform 
conditions by detecting the quasi-maximum points based on the current variation that can detect the hidden 
point. Besides the forward scanning, the proposed algorithm also tracks the GMPP using reverse scanning 
with termination criterion, reducing the search space dramatically. Comparisons with other algorithms 
developed for MPPT under PSCs are considered clarifying and showing the effectiveness of the proposed 
method. Finally, the proposed MPPT algorithm is implemented in real-time using national instruments (NI) 
CompactRIO in field-programmable gate array (FPGA) mode to demonstrate the viability of the proposed 
work. 

This paper is organized as: the effect of PSC on the PV-based power generation systems is revisited, 
and the definition of the hidden point is given in section 2. Section 3 introduces the proposed GMPPT 
algorithm. Numerical comparisons with other algorithms developed for MPPT under PSCs are provided in 
section 4. Experimental validation is discussed in section 5. Finally, conclusions are given in section 6. 


2. PV BASED POWER GENERATION SYSTEMS 
2.1. PV based power generation system characteristic under partial shading 

Several PV modulus configurations are adopted in practice to achieve the desired PV based power 
generation sizing. Without loss of generality, the PV generation system shown in Figure | is taken as an 
example to highlight the partial shading effect and define the hidden points. Because of different obstacles 
like shadows, moving clouds, and snow, the solar panels in the PV array are not facing the same weather. 
The non-uniform weather causes partial shading effects in the PV modules. The current produced by the 
shaded PV module is less than the other shine modules. This may damage the shaded module or reduce the 
generated power from the PV array. A bypass diode is connected to each module to protect the shaded 
module and minimize the power losses, as shown in Figure 1. The bypass diodes short-circuit the shaded 
modules to prevent them from operating at negative voltage [24]. However, multiple peaks exist in the P-V 
curve because of the activation of the bypass diodes Figure 2, assuming two different weather. The existence 
of multiple peaks puts stress on the adopted MPPT algorithm and may lead to tracking failure. 


2.2. The effect of the hidden points on MPPT 

Most of the GMPPT algorithms proposed in the literature use a power variation threshold (i.e., 
|AP|=5%) to detect the weather change [25]. However, monitoring the power variation alone, it is difficult to 
notice the change in weather. Figure 3 shows a PV generation system under different weathers, in which 
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Figure 3(a) shows the PV generation system under uniform weather, while Figures 3(b) and 3(c) under non- 
uniform weather. Figure 3(d) shows the P-V curves under uniform and non-uniform weather. Let the GMPP 
be located when the PV generation system faces the first partial shading condition (PSC;). The weather 
changed after some time and the PV generation system met another partial shading condition (PSC2). While 
the irradiance shifts between PSC; and PSC2, the old GMPP was preserved as a new LMPP with no power 
variation. Hence, any GMPPT algorithm relying on power variation cannot figure out this transition in 
weather. In this paper, this point is called the partial shading hidden point (PSHP). The existence of the PSHP 
reduces the maximum generated power; hence, it must be avoided. The proposed GMPPT algorithms can 
figure out the presence of PSHPs and avoid them without additional cost or burden. 
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Figure 2. The P—V and I—V curves for a shaded PV array under non-uniform weather 
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Figure 3. Series-parallel configurations PV generation system under different weathe (a) uniform irradiance 
condition (UIC), (b) partial shading condition 1 (PCS 1), (c) partial shading condition 2 (PCS2), and (d) P-V 
curves under different weather 
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3. THE PROPOSED GMPPT METHOD 

Three simple stages form the proposed algorithm: the scanning stage in which all the quasi- 
maximum power points are located and stored. The maximum of these points is passed to the second stage to 
find the exact GMPPT employing P&O. A hidden point stage is adopted to bypass the PSHP in the P-V 
curve if it exists. Moreover, the reverse scanning procedure is integrated with the proposed work to minimize 
the search space. Figure 4 shows the flowchart of the proposed GMPPT algorithm, and the detailed 
description of each stage is discussed in the following sections. 
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Figure 4. The flowchart of the proposed MPPT algorithm 


3.1. Scanning the P-V curve stage 
3.1.1. Forward scanning 

Scanning the PV curve can be either forward or reverse Figure 5. In forwarding direction scanning 
(Figure 5(a)), the reference voltage is updated from a lower voltage (Vmin) to an upper voltage (Vnax). 


Vres (kTs ) = Vref ((k — 1)Ts) + AV (1) 


where k is an integer number, T, is the sampling time, kT, is the current sample, (k — 1)T, is the previous 
sample, and AV is the voltage step. The DC-DC converter duty cycle is regulated to follow the reference 
voltage of the PV generation system. At each sample, the current is sensed, and the current variation is 
calculated using: 


I(KTs)—I((K-1)Ts) 


AI = 1((K-1)Ts) 


(2) 
The first instant at which AI > 5% is detected is the nearest to the actual peak, as shown in Figure 6. 
Therefore, an MPP is detected; this point is saved as a quasi-local maximum power point (QLMPP). The 
switch in the flowchart of Figure 5(a) is adopted for this purpose. According to Figure 5, updating the voltage 
reference is ended if the saved QLMPPs is equal to the number of series modules or the reference voltage hits 
the maximum limit. Normally, Vmin = Vino"? and Vnax = 0.9Voe™ [23]. 


3.1.2. Reverse scanning 

In the reverse scanning (Figure 5(b)), the current at Vmin is sensed and registered as Imax. Then, the 
P-V curve scanning is done from Vnay toward Vmin updating Vep as in (3) while monitoring the PV current 
as in (2). 


Vres (KTs ) = Vref ((k = 1)Ts) — AV (3) 


In the reverse scanning, the last occurrence of |AI|>5% is the closest to the actual peak, as 
demonstrated in Figure 6. Therefore, at this point, an MPP is detected and saved as a QLMPP. The switch in 
the flowchart of Figure 5(b) is adopted for this purpose. The scanning is ended if either the PV current holds 


Iny 2 0.95Imax Or the minimum voltage is reached. 
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Figure 5. Scanning the P—V curve stage adopting (a) forward scanning (b) reverse scanning 


3.2. GMPPT and partial shading detection 
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Figure 6. P-V curve forward and reverse scanning 
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In this stage, the P&O will take the maximum of the stored QLMPP to track the actual GMPP 
(Figure 7(a)). During this stage, the power variation is found using: 


AP = 


P((k-1)Ts) 


P(kTs)—P((k-1)Ts) 


(3) 


If AP > 5%, then the weather changes, and the algorithm will activate stage 3.1, scanning the P—V curve. If 
there is no variation in power for a threshold time, the algorithm will call stage 3.3, avoiding hidden points 


(Figure 7(b)). 
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Figure 7. GMPPT and PS detection (a) tracking stage (b) avoiding the hidden point stage 


3.3. Avoiding hidden point 

To prevent the PSHP, a timer is activated once the P&O tracks the actual GMPP. If a threshold time 
is over without detecting PSC observing the variation in power, the algorithm will check the saved QLMPPs, 
as shown in the flow chart of Figure 7(b). If there is any variation in the stored QLMPP power, the algorithm 
will detect a PSC and activate stage 3.1. If no change in power is noticed, the timer is initialized, and the 
algorithm will call stage 3.2, tracking and monitoring the power variation (Figure 7(a)). 


4. NUMERICAL SIMULATION 

The proposed GMPPT method in this work is independent of the adopted PV generation system 
configuration. Thus, the series-parallel configuration shown in Figure 3 is nominated without loss of 
generality. The simulated sub-modules and system parameters are given in Table 1. Figure 8 shows the 
implementation of the proposed GMPPT algorithm in which a buck-boost DC-DC converter interfaces the 


PV generation system to a load. 
The Proposed Vref 
GMPPT 
Algorithm o] 


Table 1. The sub-module and simulation 


parameters 

Parameters Symbole Value I 
Output power Pax 35 W 
Open circuit voltage Vinod 11.1 V 
Short circuit current jae 415A 
Voltage at Prax Vinpot 9V 
Current at Prax proa 3.89 A 
Buck-Boost inductance L limH EOE ~ Buck = Boost converter 
Input side capacitance Cin 470 uF 
Output side capacitance Cüt 100 uF ? : : 
Switching Frequency f, 20 kHz Figure 8. The implementation of the proposed GMPPT 
Resistive load R, 100 algorithm 


4.1. Hidden point simulation 

This simulation verifies the capability of the proposed GMPPT algorithm to detect and avoid the 
PSHP. The forward scanning-based algorithm (Figure 5(a)) is nominated without loss of generality. The 
simulation will mimic the scenarios shown in Figure 3 and uses the system shown in Figure 8. Three kinds of 
weather are taking place: a UIC for one second, PSC1 for another one second, and PSC2 for the remaining 
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time of the simulation. The proposed algorithm starts with the scanning stage and finds one peak (P = 
850 W) at (V = 80 V). The P&O finds the actual GMPP (P = 853.3 W) at (V = 81.7 V) while monitoring 
the power. The hidden point stage of the proposed method is activated once the timer hits the timer threshold. 

Under UIC, no change in power has been detected; hence, the algorithm returns to stage 3.2. A PSC 
is detected at 0.9 sec.; hence the algorithm initiates stage 3.1, locating two QLMPP. The maximum among 
them is passed to stage 3.2, and the P&O finds the GMPP. At 1.4 sec., the threshold time is reached, and 
stage 3.3 is activated. The algorithm examines the saved QLMPP; no power variation is noticed, indicating 
no weather change. Therefore, the timer is initialized, and stage 3.2 is triggered. A change in the weather 
occurs at t =1.9 sec., with an LMPP equal to the previous GMPP. Hence, any algorithm that relies on the 
power variation will fail. However, the proposed algorithm can detect this event by checking the saved 
QLMPPs power. A change in power is noticed. Thus, stage 3.1 is called. Figure 9 shows the results that 
clarify and demonstrate the proposed method’s capability to overcome the PSHPs. 
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Figure 9. The hidden point simulations 


4.2. Comparisons with other methods 

Comparisons with the works given in [22], [23] are carried out to demonstrate the effectiveness of 
the proposed GMPPT algorithm. Figure 10 shows the P-V curve under different weather conditions. At first, 
the UIC1 occurs between te[0,0.59) sec. At t=0.59 sec, PSC1 occurs until t=1.47 sec, followed by PSC2 
until 2.2 sec, and UIC2 until =3.5 sec. Adopting the forward scanning, all the voltage range between Vmin 
and Vnax is scanned, while it is unnecessary to examine all the voltage range in the reverse scanning method. 
Figure 11 shows the browsed voltage range adopting the reverse scanning under each condition. Figure 12 
shows the comparison between the proposed and other algorithms, in which Figure 12(a) depicts the results 
of the proposed method adopting the forward scanning. In contrast, Figure 12(b) shows the results of the 
proposed method assuming reverse scanning. Figures 12(c) and 12(d) show the results adopting the 
approximate I-V and the PV current variation-based methods, respectively. Figure 13 shows a zoom-in of 
comparing the proposed and other algorithms. Figures 12 and 13 show that although all utilized methods can 
track the GMPP, the proposed method has the shortest tracking time. Hence, the higher average power is 


generated adopting the proposed method. A summary of the results is given in Table 2. 
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Figure 10. P-V curves under different PSC 
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Figure 11. The browsed voltage range adopting the reverse scanning 


Array Voltage(V) 


Array power(w) 
o 8 88 2 2 


` 3 0 0.5 1 1.5 2 2.5 3 
Time(s) Time(s) 


Proposed method adopting reverse scanning 


Array Voltage(V) 
D a 
è 8 


iS) 
© 


Array power(w) 
2 £228 
>o 8 8 83 8 8 


0 0.5 1 1.5 2 2.5 3 0 0.5 1 1.5 2 2.5 3 
Time(s) Time(s) 


Method in [23] 


e] 


Array power(w) 
[> pS a æ z 
o § 8 8 8 8 


D 


Array Voltage(V) 
n 


w 


0 0.5 1 1.5 2 25 3 1.5 2 2.5 3 
Time(s) Time(s) 


o 
fsd 
D 


Method in [22] 


g 


Array Voltage(V) 
S 


Array power(w) 


w 


0 0.5 1 1.5 2 2.5 


Time(s) Time(s) 


(d) 


Figure 12. Comparison between the proposed and other algorithms (a) the proposed method adopting the 
forward scanning, (b) the proposed method adopting the reverse scanning, (c) approximate I-V based 
method, and (d) PV current variation-based method 
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Figure 13. Zoom-in of the comparison between the proposed and other algorithms 


Table 2. Comparison between the proposed and other algorithms developed for MPPT under PSCs 


Method Total Tracking time Average power Storing LMPPs Calling the P&O 
The proposed work integrating forward scanning 140 ms 571.64 watt Yes Once 
The proposed work integrating reverse scanning 90 ms 573.33 watt Yes Once 
The method in [23] 200 ms 570.84 watt No Once 
The method in [22] 294 ms 567.77 watt Yes Number of LMPPs 


5. EXPERIMENTAL VALIDATION 

The LabVIEW CompactRIO 9063 is used to validate the proposed method experimentally, as shown 
in Figure 14. The current is sensed using shunt resistance and interfaced to NI9225 analog to digital converter. 
The PV voltage is interfaced directly to the NI9225. Three PV modules (PS M36s-100W) manufactured by 
Philadelphia Soler are connected in series. The used PV module has the data as listed in Table 3. 
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Figure 14. Experimental setup test-bed 


Table 3. Data of PV module PS M36s-100W 


Parameters Value 
Papp 100 +3% W 
Voc 22.64 V 
Toe 5.72A 
Vinpp 18.39 V 
Impp 5.44 A 
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Four scenarios are considered to validate the proposed work; the data for each scenario has been 
collected and listed in Table 4. The first scenario has no partial shadow; in this case, the proposed algorithm 
detects only one peak at 50.17 V, and the P&O locates the exact GMPP at 47.86 V. One of the PV modules 
faces a partial shadow in the second scenario. The scanning stage of the proposed method locates two 
QLMPPSs at 41.07 and 55.25 V, respectively. The maximum at 41.07 V is given to the tracking stage, and the 
P&O finds the exact GMPP at 39.09 V. 

The third test is for two partially shaded modules. In this case, three peaks are detected; the 
maximum among them is given to the P&O to find the exact GMPP. Finally, all of the modules are facing 
different irradiances. In this case, three peaks are detected; the maximum among them is given to the P&O to 
find the exact GMPP. The results listed in Table 4 validate the applicability of the proposed MPPT algorithm 
in locating the GMPP under partial shading conditions. 


Table 4. Experimental test data under different irradiances 


Partial Shading condition LMPP, LMPP; LMPP; GMPP 
PV, PV, PV; Vi P; V2 Pı V3 P3 Vp&o Pp&o 
— — — 50.17 182.29 47.86 191.23 
— PS — 41.07 150.51 55.25 115.37 — — 39.09 156.7 
PS PS — 30.69 119.69 44.89 108.67 55.23 114.12 30.25 121.08 
PS PS PS 21.58 87.07 44.92 108.22 55.2 113.44 54.49 115.4 


6. CONCLUSION 

An algorithm to extract the maximum power from the PV-based generation systems under non- 
uniform weather, with the ability to detect and bypass the hidden point in the partially shaded P—V curve, is 
presented in this paper. PSHP is caused by a transition between GMPP to an LMPP once the partial shading 
patterns have been changed, which is hard to be observed monitoring the power difference alone. Three 
simple stages form the proposed algorithm: the scanning stage in which all the quasi-maximum power points 
are located and stored. The maximum of these points is passed to the second stage to find the exact GMPPT 
employing P&O. A hidden point stage is adopted to bypass the PSHP in the P—V curve if it exists. The 
reverse scanning procedure is integrated with the proposed work to minimize the search space. The 
performance indices of the proposed GMPPT algorithm, compared with other algorithms developed for 
MPPT under PSCs, are better in terms of the generated power and the time to track the actual global peak. 
The proposed MPPT algorithm is two times faster than the compared methods, generating almost 2% extra 
power. Besides the simulation verification, the proposed MPPT algorithm is implemented in real-time using 
NI Compact-RIO in FPGA mode to demonstrate the viability of the proposed work. 
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